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Q> ' Abstract. The power spectrum of mass density fluctuations is evalu- 

CN \ ated from the Mark III and the SFI catalogs of peculiar velocities by 

a maximum likelihood analysis, using parametric models for the power 
spectrum and for the errors. The applications to the two different data 
' sets, using generalized CDM models with and without COBE normal- 

OO . ization, give consistent results. The general result is a relatively high 

^ ' amplitude of the power spectrum, with a%VLr^'^ = 0.8 ± 0.2 at 90% con- 

fidence. Casting the results in the 0^ — ^^A plane, yields complementary 
I constraints to those of the high-redshift supernovae, together favoring 

. a nearly flat, unbound and accelerating universe with comparable con- 

tributions from and Oa- Further implications on the cosmological 
Q^. parameters, arising from a joint analysis of the velocities together with 

' \ small-scale CMB anisotropics and the high-redshift supernovae, are also 

^ ■ briefly described. 

to 
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V~j ■ 1. Introduction 

■ In the standard picture of cosmology, structure evolved from small density fluc- 

tuations that grew by gravitational instability. These initial fluctuations are 
assumed to have a Gaussian distribution characterized by the power spectrum 
(PS). On large scales, the fluctuations are linear even at late times and still 
governed by the initial PS. The PS is thus a useful statistic for large-scale struc- 
ture, providing constraints on cosmology and theories of structure formation. 
It is advantageous to estimate it from velocities, as these are directly related 
to the mass density fluctuations, and are effected by large scales and thus are 
approximately still linear. 

In this work, we develop and apply a likelihood analysis for estimating the 
mass PS from peculiar velocities. This method uses the "raw" peculiar veloc- 
ities without additional processing, and thus utilizes much of the information 
content of the data. It also takes into account properly the measurement errors 
and the finite discrete sampling. The simplifying assumptions made are that 
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the velocities follow a Gaussian distribution and that their correlations can be 
derived from the density PS using linear theory. 

We use the two comprehensive available peculiar velocity data sets, the 
Mark III catalog (Willick ct al. 1997) and the SFI catalog (Haynes et al. 1999a,b; 
Giovanclli, these proceedings). Mark III samples ~ 3000 galaxies within a dis- 
tance of ~ 70/i^^Mpc around us, and SFI consists of ~ 1300 spiral galaxies with 
a more uniform spatial coverage in a similar volume. The typical relative dis- 
tance errors of individual galaxies are 15 — 20%, and both catalogs are carefully 
corrected for systematic errors. It is interesting to compare the results of the 
two catalogs, especially in view of apparent discrepancies in the appearance of 
the velocity fields (e.g., da Costa et al. 1996, 1998). We explore the cosmological 
implications of the velocities on their own and in conjunction with constraints 
derived from other types of data. 

2. Method 

Given a data set d, the goal is to estimate the most likely model m. Invok- 
ing a Bayesian approach (and assuming a uniform prior), this can be turned 
to maximizing the likelihood function L = V{d\m), the probability of the 
data given the model, as a function of the model parameters. Under the 

assumption that both the underlying velocities and the observational errors 
are Gaussian random fields, the likelihood function can be written as £ = 
[(27r)^ det(i?)]-i/2 g^p (-1 Y,fj diRr^^dj) , where {dJili is the set of observed 

peculiar velocities and R is their correlation matrix. R involves the theoretical 
correlation, calculated in linear theory for each assumed cosmological model, 
and the estimated covariance of the errors. 

The likelihood analysis is performed by choosing some parametric functional 
form for the PS. Going over the parameter space and calculating the likelihoods 
for the different models, one finds the PS parameters for which the maximum 
likelihood is obtained. Note that this method, based on velocities, in fact mea- 
sures P{k)^jn^''^ and not the PS by itself. Confidence levels are estimated by 
approximating — 21n>C as a distribution with respect to the model parame- 
ters. We have extensively tested the reliability of the method on realistic mock 
catalogs, designed to mimic in detail the real data. 

Our main analysis is done with a suite of generalized CDM models, normal- 
ized by the COBE 4-year data. These include open models, flat models with a 
cosmological constant and tilted models with and without a tensor component, 
where the free parameters are the mass-density parameter Qrn, the Hubble pa- 
rameter h and the power index n. We also use CDM-like models where the 
amplitude was allowed to vary. 

A common problem in PS estimations is that the recovered PS is sensitive 
to the assumed observational errors, that contribute as well to the correlation 
matrix R. To alleviate this problem, we extend the method so that also the 
magnitude of these errors is determined by the likelihood analysis. This is done 
by adding free parameters that govern global changes of the assumed errors, in 
addition to modeling the PS, and provides a consistency check of the magnitude 
of the errors. We find, for both catalogs, a good agreement with the original 
error estimates, which allows for a more reliable recovery of the PS. 
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Figure 1. Likelihood analysis results for the COBE-normalized flat 
ACDM model with h = 0.6. Shown are ln£ contours in the Qm — n 
plane for SFI (top left panel; cf. Frcudling et al. 1999) and for Mark III 
(top right; cf. Zaroubi et al. 1997). The best-fit parameters for SFI 
and Mark III are marked on both by 'S' and 'M', respectively. The 
lower panel shows the corresponding maximum-likelihood PS for SFI 
(solid line) and for Mark III (dashed). The shaded region is the SFI 
90% confidence region. The three solid dots mark the PS calculated 
from Mark III by Kolatt and Dekel (1997), together with their quoted 
la error-bar. (in the plot = ^m) 
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3. Results 



Figure |^ shows, as a representative example, the results for the COBE-normalized 
flat ACDM family of models, with a tensor component in the initial fluctua- 
tions, when setting h = 0.6 and varying Q^n and n. Shown are InC contours 
for the SFI catalog and for Mark III. As can be seen from the elongated con- 
tours, what is determined well is not a specific point but a high likelihood ridge, 
constraining a degenerate combination of the parameters roughly of the form 
J^„n3-7 = 0.6 ±0.1, in this case. The corresponding best-fit PS for the two cata- 
logs is presented as well, with the shaded region illustrating the 90% confidence 
region obtained from the SFI high-likelihood ridge. 

These results are typical for all other PS models we tried. For each catalog, 
the different models (including the ones with free normalization) yield similar 
best-fit PS, falling well within each others formal uncertainties and agreeing 
especially well on intermediate scales {k ~ 0.1/iMpc~^). The similarity of the 
PS obtained from SFI with that of Mark III, which is apparent in the figure, is 
also illustrative of the other models. This indicates that the peculiar velocities of 
the two catalogs, with their respective error estimates, are consistent with arising 
from the same underlying mass density PS. This does not preclude possible 
differences that are not picked up by this statistic, but can be viewed as another 
indication of the robustness of the results. Note also the agreement with an 
independent measure of the PS from the Mark III catalog, using the smoothed 
density field recovered by POTENT (the three dots; Kolatt & Dekel 1997). 

More details on the method and results can be found in Zaroubi et al. (1997) 
and Freudling et al. (1999), regarding the applications to Mark III and to SFI, 
respectively. The robust result, for both catalogs and all models, is a relatively 
high PS, with e.g. P{k = 0.1 /iMpc-^)J7m^-^ = (4.5 ± 2.0) x 10^ {h'^Mpc)^. An 
extrapolation to smaller scales using the different CDM models gives asQm^'^ = 
0.8 ± 0.2. The high-likehhood ridge is a feature of all COBE-normahzed CDM 
models, corresponding to a general constraints on the combination of cosmologi- 
cal parameters of the sort ilm h^o^ = 0.6±0.2, where ^ = 1.3 and v = 3.7, 2.0 
for flat ACDM models with and without tensor fluctuations respectively. For 
open CDM, without tensor fluctuations, the powers are ^ = 0.9 and v = 1.4. 
These error-bars are crude, reflecting the 90% formal likelihood uncertainty for 
each model, and the variance among different models and catalogs. Note that 
our results are suggestive of somewhat higher values of a^Vtm^'^ and than 
those implied by some other methods (such as from the clusters abundance and 
the different /3 measures). 

4. — r^A Constraints 

We have recently extended the analysis of COBE-normalized CDM models to 
models with general values of Vim and Q.\ (Zehavi & Dekel 1999). Although 
the Oa dependence comes in only indirectly through the COBE normalization, 
such results are particularly interesting as they can be combined with other 
constraints in the 0^ — ^^A plane. 

Figure ^ illustrates constraints in the ^.^ — ^k plane, showing the confidence 
contours obtained by the Supernova Cosmology Project (Perlmutter et al. 1999; 
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Figure 2. Constraints in the — plane arising from the high- 
z SN (the blue inclined contours; Perlmutter et al. 1999) and from 
a likelihood analysis of the SFI peculiar velocities (the red roughly 
vertical contours). The 68, 90 and 99% confidence regions are shown 
for both. The peculiar velocity contours are for n = 1, h = 0.65, 
with the shifted red dashed lines showing the effect of changing the 
values of these parameters. The corresponding joint confidence regions 
of the velocities and the SN are shown as the black overlapping ellipses 
(Zehavi & Dekel 1999). 
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consistent with the findings of the High-z Supernova Search Team, Riess et 
al. 1998) as well as the contours from our SFI hkehhood analysis, performed 
for fixed values of n = 1 and h = 0.65. The velocity analysis constrains an 
elongated ridge in this plane of a nearly fixed 0,^ and varying The analogous 
constraints from the Mark III catalog are quite similar, except that the (fairly 
uncertain) upper bounds on $7a are slightly tighter. A change in the values of n 
and h essentially shifts the ridge toward a higher or lower 0.^, for smaller and 
larger values of these parameters, respectively. (This is another manifestation of 
the degeneracy between these parameters mentioned earlier). Their acceptable 
range is therefore needed to be determined a-priori from other constraints. For 
a reasonable range of these parameters, the effect is denoted in the plot by the 
shifted dashed lines. 

The corresponding joint contours of the velocities and the SN, obtained 
by multiplying the likelihoods, are also shown on the plot, and here as well one 
should consider the area bounded by the dashed lines. Taking into consideration 
concurrently these two independent sets of constraints thus implies a consider- 
able contribution from both and Oa- Specifically, models with small and 
small r^A, that are still allowed by the SN constraints alone, are disfavored when 
the constraints from SNe and peculiar velocities are considered jointly; together, 
they make a stronger case for an unbound accelerating universe with a positive 
cosmological constant. 

5. Further Analysis 

We are currently in the final stages of performing a joint analysis of peculiar 
velocities together with constraints obtained from CMB anisotropics and from 
the high-z SN (Bridle et al. 1999). The distinct types of data complement one 
another, each constraining different combinations of the cosmological parame- 
ters, and together can potentially set tight constraints. Figure ^ illustrates the 
constraints obtained for each of these data sets in the as — h — space, for 
the scale-invariant flat ACDM family of models. The results shown here are for 
the SFI catalog. Since the COBE constraint is now included in the CMB data, 
the velocity analysis is performed free of the COBE normalization. The high-z 
SN constraints appear here as bounds on (as we assume flat ACDM). The 
value of h is left free to be determined by the likelihood analysis. 

The figure demonstrates that these three data sets provide roughly orthog- 
onal constraints, with a significant overlap beyond the 2cr level, that allows a 
meaningful joint parameter estimation. The best-fit values obtained from the 
joint analysis of the three data sets are = 0.52, h = 0.57 and erg = 1.1 
(corresponding to ag^m^'^ = 0.74). The 95% confidence regions on the individ- 
ual parameters, obtained by marginalizing over the other two parameters, are 
0.27 <nm< 0.54, 0.54 <h< 0.85 and 1.02 < erg < 1.67. See also Lahav, these 
proceedings, for a presentation of the results in the ag^m^'^ — ^mh plane. 

Additional work in progress includes an application of the likelihood method 
to other peculiar- velocity data, such as velocities of galaxy clusters (e.g., the 
SMAC sample; Hudson et al. 1999, Smith, these proceedings) or the local SN 
velocities, both probing larger scales with relatively high accuracy per object. 
With regard to the cluster velocities, we are attempting to properly take into 
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Figure 3. The 95% confidence regions obtained from the SFI veloci- 
ties, the CMB anisotropies and the high-z SN, in the erg — /i — space, 
for scale- invariant flat ACDM models (Bridle et al. 1999). 
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account in the analysis the fact that they are sampled at high density regions 
(Zehavi et al., in prep.). 

Finally, another project underway is an attempt to obtain model-independent 
band-power estimates of the PS from peculiar velocities, using an iterative 
quadratic estimator scheme, which greatly improves the computational effort. 
(Such an approach is commonly applied to CMB measurements, e.g Bond, Jaffe, 
&; Knox 1998). This allows to relax the a-priori assumption of the PS form, and 
would illuminate the actual constraints on the different scales and their cross- 
correlations (Zehavi &: Knox 1999). 
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